Compact optical circulator based on a uniformly magnetized ring cavity 
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We propose a new class of compact integrated optical circulators providing a large isolation level 
while maintaining a straightforward technological feasibility. Their layout is based on a nonrecip- 
rocal radial Bragg cavity composed of concentric magneto-optical rings. The circulator ports are 
standard rib waveguides, butt-coupled to the cavity by cutting through its outer rings. The device 
is specifically designed for operation in a uniform external magnetic field. Using a coupled-mode 
description of the complete cavity/waveguide-port system, we explore the rich behaviour of cavity 
circulators in presence of varying levels of direct port-to-port coupling. We demonstrate numerically 
a strongly miniaturized two-dimensional cavity circulator, with a total footprint of less than (fOA) 2 , 
achieving a 20-dB isolation level at telecom frequencies over a bandwidth of f 30 GHz. The device 
is found to be very tolerant with respect to fabrication imperfections. We finish with an outlook on 
three-dimensional versions of the proposed nonreciprocal cavities. 
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I. INTRODUCTION 



The need for an integrated and miniaturized version of 
an optical isolator, or more generally an optical circula- 
tor, is making itself increasingly felt. The drive towards 
ever higher degrees of all-optical on-chip integration is 
often hindered by the absence of an element that induces 
a one-way sense in the path of the signal in such an in- 
tegrated circuit. Without it, long path interferences can 
lead to large changes somewhere in the circuit due to 
small amplitude oscillations at a remote point. Com- 
mercially available isolators are bulk free-space devices 
based on 45° nonreciprocal magneto-optical (MO) Fara- 
day polarization rotators in combination with polariz- 
ers placed at their entrance and exit. The most com- 
monly used MO materials are magnetic garnet oxides 
such as Ce-substituted Ce x Y^-xFe^Oi2 (Ce:YIG), com- 
bining optical transparency and strong MO properties at 
telecom frequencies. Realisation of an integrated isola- 
tor based on the Faraday-effect is very difficult because 
of the inevitable geometric birefringence of planar inte- 
grated waveguide circuits. 

Nowadays, research efforts focus on integrated isola- 
tor concepts that do not depend on polarization con- 
version [1]. This can be achieved by properly orienting 
the magnetization in the MO material so that no cou- 
pling occurs between the quasi- TE and quasi-TM wave- 
guide modes. Phase velocity and field profiles of the 
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waveguide modes then become different for forward and 
backward propagation, while their polarization state is 
unchanged. This has been exploited to propose vari- 
ous concepts based on nonreciprocal interference [2, 3], 
multimode imaging [4], microring and -disk resonators 
[5, 6], etc. Nevertheless, experimental demonstration of 
garnet-based isolators suffers from the limited MO gy- 
rotropy in the optical and near-infrared regime. This 
leads to high device lengths (typically of the order of 
1 mm) , which goes hand in hand with serious technolog- 
ical challenges, for instance to maintain magnetic uni- 
formity. As a result, during the past decade the use of 
resonant photonic-crystal (PhC) layouts has attracted a 
lot of interest. Magnetophotonic crystals allow to arti- 
ficially enhance the intrinsic strength of the basic MO 
effects [7, 8]; in addition, removal of time-reversal sym- 
metry in periodic structures leads to the appearance of 
novel phenomena, such as frozen light [9] , photonic chiral 
edge states [10] and one-way band gaps [11, 12]. Over the 
last few years, an increasing number of promising minia- 
turized isolator (and circulator) designs using PhC effects 
have been reported [11, 13-16]. 

In this article, we will focus on a particular class of 
miniaturized integrated circulator designs based on the 
resonant enhancement of the light-MO-material interac- 
tion in a nonreciprocal cavity. Uniformly magnetized 
resonant ferrite cavities have been used for decades in 
microwave circulators [17]. However, optical circula- 
tors made by a simple geometric rescaling of existing 
microwave devices would have prohibitively low opera- 
tion bandwidths. This is because at optical frequencies 
nonreciprocal effects arc induced by the gyroelectric off- 
diagonal elements of the permittivity tensor, which are 
typically one or two orders of magnitude smaller than the 
analogous gyromagnetic off-diagonal elements of the per- 
meability tensor in the microwave regime [18]. In 2005, 
Wang and Fan reported a solution for a cavity circulator 



operating at optical frequencies [13]. The proposed de- 
vice is composed of a 2D PhC cavity etched in bismuth 
iron garnet (BIG) — a transparent magnetic oxide with 
record MO properties [19]. The cavity is evanescently 
coupled to three symmetrically placed PhC waveguides. 
It achieves an infinitely strong circulation at the reso- 
nant wavelength of the cavity and has a footprint of just 
a few square wavelengths. However, under uniform mag- 
netization the spectral bandwidth of this circulator be- 
comes negligibly small, reducing ultimately the applica- 
bility of the device. A reasonable bandwidth of the order 
of 50 GHz can only be achieved by imposing a very spe- 
cific domain structure of antiparallel magnetic domains. 
Achieving and maintaining this magnetic domain struc- 
ture within an area of a few jim 2 is unfeasible. 

In order to remedy the unfeasibility of this concept we 
reported earlier an original design approach for a PhC 
cavity that provides simultaneously high circulation lev- 
els (> 30 dB), good spectral bandwidth (~80GHz), and 
operates in a uniformly magnetized MO material [14]. 
The new concept is based on an approximately axisym- 
metric arrangement of the holes making up the PhC cav- 
ity. In this work we propose to eliminate entirely the use 
of the PhC layout, forming instead a fully axisymmet- 
ric cavity composed of concentric MO rings and butt- 
coupling it to standard rib waveguides. The general ge- 
ometry of the proposed new cavity circulator is shown 
in Fig. 1. If this device shows similar performance as 
the cavity circulators proposed in Refs. 13, 14 and 5, its 
markedly simpler layout would make it largely preferable 
over the existing circulator schemes. There are neverthe- 
less important issues to be tackled. The circular Bragg 
grating formed by the concentric full and split rings can in 
theory provide sufficient confinement [20, 21]. However, 
the quality factor will be lower than that of a PhC cavity. 
Moreover, this ring circulator is not embedded in a PhC 
crystal and thus does not operate in a bandgap. As a 
result, its operation can be perturbed by power loss to 
free space and direct waveguide-to-waveguide coupling. 
This latter mechanism allows the appearance of Fano res- 
onances in the transmission spectrum of the cavity. In 
this paper we will demonstrate that, even in the pres- 
ence of these effects, the circulator layout of Fig. 1 can 
be designed to achieve competitive performance. 

This paper is organized as follows. We begin by ex- 
tending, in Section II, the temporal coupled-mode model 
introduced by Wang and Fan [13] to the case of non- 
negligible direct waveguide-to- waveguide coupling and in- 
plane radiation losses. This will allow a first assessment 
of the theoretical maximally achievable device perfor- 
mance. Section III will then focus on the design and 
performance of the new type of cavity circulator. We 
will first detail how the ring cavity itself is optimally 
designed to allow operation in a uniform magnetic field 
(Section III A), and then deal with the design optimiza- 
tion of the complete circulator (Sections III B and III C), 
using the insight provided by the coupled-mode model 
of Section II. In this way it will be shown how this type 
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FIG. 1 Geometry of a circulator composed of a ring cavity 
butt-coupled to three rib- waveguides. A uniform magnetic 

field perpendicular to the xy plane causes nonreciprocal 
coupling of the even and odd degenerate if-polarized cavity 
modes. The number of split and full rings, n s resp. rif, the 
width of the slits in the split rings, d^n, and the distance 
from the waveguide ends to the centre of the cavity, p wg , arc 
all design optimization parameters. 



of circulator competes with the PhC cavity layouts of 
our earlier work [14]. Before concluding, in Section IV 
we will indicate the additional constraints placed on the 
cavity design by the presence of a super- and substrate, 
and show initial designs of properly behaving realistic 3D 
versions of circulator cavities. 



II. TEMPORAL COUPLED-MODE THEORY OF 
RING-CAVITY CIRCULATORS 

In contrast to PhC-based circulators, where the 
bandgap provides an almost perfect discoupling of the 
waveguides, in systems based on rib waveguides one can- 
not a priori exclude the possibility of direct waveguidc- 
to-waveguide coupling. In this section we will firstly ex- 
tend the abstract temporal coupled- wave description of a 
nonreciprocal resonator to include such direct transmis- 
sions (Section II A), using the work of Suh et al. [22] as 
a basis. Once the model equations have been derived, 
their formal solution, in Section II B, will allow to iden- 
tify the critical device parameters and to evaluate the 
possible detrimental influence of the direct pathway cou- 
pling on the circulator performance. As will be shown 
in Section II C, nonreciprocal Fano- type resonances ap- 
pear in the circulator transfer function when the device's 
ports are not perfectly decoupled. This does not nec- 
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essarily undermine the device's behaviour, however, but 
even allows novel functionalities. In the last subsection, 
II D, radiative cavity decay will be included in the model 
and its influence on the circulator performance will be 
quantitatively evaluated. 



A. Inclusion of direct pathway coupling 

Let us consider a circulator composed of three iden- 
tical single- mode waveguides, Wi, W2, and W3, weakly 
coupled with a resonant cavity and arranged so that the 
whole system has C% v symmetry. A particular example 
of such a system is shown in Fig. 1. In the absence of 
MO coupling, the cavity is assumed to support a pair 
of degenerate orthonormal eigenmodes of frequency co , 
belonging to the unique two-dimensional irreducible rep- 
resentation of the Czv point group. They can be classified 
as even or odd, according to their symmetry with respect 
to reflection about the axis of waveguide W\, which is 
assumed to lie along the — x axis. The circulator's op- 
eration at frequency lo is described by the coupled-mode 
equations [22] 



-vjja = — (ifi + f )a + D T s- ln , 
«out = Csi n + Da, 



(la) 
(lb) 



where the vector a = (a c ,a ) T contains the amplitudes 
of the even and odd cavity mode, the vectors s; n = 
(si.in, s 2 ,in, S3,in) T and s out = (si, ou t , S2,out , s 3 , ut ) T con- 
tain the amplitudes of the in- and outgoing waveguide 
modes, the matrix C describes the direct coupling be- 
tween waveguides, D the cavity- waveguide coupling, Q 
the cavity mode eigenfrequencies and their coupling, and 
T their decay. These matrices are subject to the funda- 
mental constraints [22] 



D^D = 2f, 

CD* = -A 



(2a) 
(2b) 



where * denotes complex conjugation and t Hcrmitian 
conjugation. Equation (2a) means in particular that all 
the energy of the decaying cavity modes is passed to the 
waveguides, i.e., the system is lossless as a whole. The 
case with radiation loss will be investigated later. To- 
gether with those following from the system's symmetry, 
the above constraints can be used to reduce the number 
of independent parameters necessary for the characteri- 
zation of the device. 

The time-reversal symmetry breaking of the device is 
assumed to be entirely described by the cross-coupling of 
the cavity resonances via the off-diagonal elements of the 
Hcrmitian matrix f2. In Section III A we show that these 
are induced by magnetization perpendicular to the cav- 
ity's plane, which couples nonreciprocally its p-polarized 
modes (with the electric field oriented in-plane). The 
coupling to the waveguide output ports is considered to 



be reciprocal, as testified by the use of the same cavity- 
waveguide coupling matrix D in both Eqs. (la) and (lb). 
Similarly, the constraint described by Eq. (2b) presup- 
poses time-reversal symmetry of the cavity-waveguide 
coupling [22]. This ansatz is justified since the main MO 
interaction in the device is the nonreciprocal coupling of 
the cavity resonances among themselves. We shall now 
consider each of the matrices occurring in Eq. (1) in turn. 

Owing to the three-fold rotational symmetry of the cir- 
culator, the C matrix must have the form 



C = 



r t t 
t r t 
t t r 



(3) 



where t = \t\ e 1T and r = |r| c 1 ( 7r+T+A ) arc the transmis- 
sion and reflection coefficients of the waveguide modes 
in the absence of the cavity. For future convenience we 
include an explicit term tt in the phase factor of r. All 
nonreciprocity is supposed to be contained in the cavity 
mode coupling, hence the symmetry of C. If the system 
is lossless, as we assume in the present subsection, we can 
use the condition of unitarity of C (or |r| 2 +2|t| 2 = 1 and 
\t\ — 2\r\ cos A = 0) to express t and r as 



t = 



2cosAe iT 
Vl + 8cos 2 A' 



pi(r+A) 



Vl + 8cos 2 A 



(4) 



with the additional condition cos A > 0. Without loss 
of generality we can therefore restrict A to the interval 
[— j/7r, jj7r]. Note that the case of no direct coupling, con- 
sidered in Ref. 23, corresponds to A = ±^tt (in other 
words, the direct waveguide transmission and reflection 
coefficients being in quadrature), while maximum direct 
coupling occurs when A = (r and t in antiphase). 
The D matrix has the general form 



D 



di e di 
d2c d 2o 
die d 3o 



(5) 



where d im (i = 1, 2, 3; m — c, o) describes the coupling of 
mth cavity mode with ith waveguide. These coupling pa- 
rameters are proportional to the values of the electromag- 
netic field of the modes along the waveguide axes. Ow- 
ing to the assumed symmetry of the mode fields and the 
symmetrical arrangement of the three waveguide ports, 
these 6 complex parameters can be expressed in terms of 
a single complex coupling constant d = \d\ e lS : 



D = d 



1 ' 

_ 1 V3 

2 2 

_ 1 _V3 

2 2 ■ 



= dD. 



(6) 



Substituting this formula into Eq. (2a), we obtain 



1 
1 



, where 7 = ||d| 2 . 



(7) 
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The second constraint, Eq. (2b), yields 
(t-r)d*=d, or d 2 = (t-r)\d\ 2 



7 (i-r). (8) 



The latter form is important, since we shall see further 
down that the behaviour of the device depends solely 
on d 2 . This expression therefore shows that the phase 8 
of the coupling constant d plays no role. From Eq. (4) 
we have 



t-r = 



2 cos A + e iA 
Vl + 8cos 2 A 



3 cos A + i sin A 
13 cos A + isin A| 



(9) 



gi[r+arg(3 cos A+i sin A)] 



The coupling between the cavity modes is described 
by the off-diagonal elements in Cl. Since we presuppose 
nonreciprocal and lossless coupling, the Cl matrix must 
take the form [23] 



n = 



w v 
-V uj 



(10) 



where V = iV is purely imaginary. The real quan- 
tity V will henceforth be called the (MO) mode coupling 
strength. Its form will be detailed in Section III A. The 
eigenvalues of the above matrix, uj± = uj a ± V, are the 
frequencies of the eigenmodes of the cavity in isolation 
(uncoupled to waveguides). The corresponding eigen- 
vectors are obviously (l,±i) T , i.e. counterclockwise- and 
clockwise-rotating combinations of the even and odd cav- 
ity modes. The frequency splitting Alu = \lo + — cj_| = 
2\V\ is proportional to the mode coupling strength. 



B. Solution of the coupled- mode model 

We have now collected enough information to solve the 
coupled-mode equations (1) and describe the device be- 
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haviour in terms of a minimal number of parameters. 
These are (1) the frequency cj of the degenerate cav- 
ity resonances, (2) the MO mode coupling strength V, 
(3) the decay rate 7 [related via Eq. (7) to the cavity- 
port coupling strength] , (4) the phase difference A + n 
between the scattering coefficients r and t, and finally 
(5) the phase r of the direct pathway waveguide trans- 
mission coefficient t. 

The formal solution of the coupled-mode equations (1) 
reads 



Sout = {C + D[i((l - tui) + f]- 1 ^}^ 



(11) 



It is already apparent that device operation only depends 
quadratically on the coupling coefficient d. Taking W\ to 
be the input waveguide by setting s; n = (1,0, 0) T , and 
noting that (D T )*i, the first column of the matrix D T , 
equals d(l,0) T , one obtains 



sout - di + d 2 {D[i(n - ui) + f]- 1 }^. (12) 

The inverse of the matrix in brackets exists if and only if 
7 7^ or u> 7^ cj ± V, and is equal to 



- i(uj — ljo) V 
— V 7 — i(w — coo) 

(13) 

We can now use Eqs. (3), (4), (6), (8), (9) and (13) to 
express the quantities occurring in Eq. (12) in terms of 
the five parameters listed in the previous paragraph. This 
leads to 



[ 7 -i( W - Wo )] 2 + ^ 2 



Sl,out 



S2,out 



S3,out 



VI + 8 cos 2 A 



VI + 8 cos 2 A I 



iA 4(3 cos A + isin A) 7 [7 — i(w — cj )] 



2e' 



VI + 8 cos 2 A I 



cos A 



3 [7- 


- i(w - 


-uo)} 2 - 


f V 2 


3 cos A + i sin A 7(7 + 


vVs 


— 'i(lo — 


wo)] 


3 [ 7 _ 


i(u> — 


uo)} 2 + 


V 2 


3 cos A + i sin A 7(7 — 


VV3 


— i(uj — 


w o)] 



[7 - i(w - c^o)] 2 + V 2 



(14a) 
(14b) 
(14c) 



The transmission amplitudes S2,out and S3 j0 ut trans- 
form into each other when the sign of the MO coupling 
strength V is inverted. This is consistent with the fact 
that V is proportional to the magnetization of the MO 
material, as will be proven in Section III A, and that, 



as an axial vector, magnetization is inverted when mir- 
rored in a plane parallel to it [18]. As a result the W3- 
output of the ring circulator in Fig. 1 must equal the 
(mirrored) WVoutput with opposite magnetization (and 
thus V), and vice versa. It is easily seen that the re- 
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flectance R = |si iOUt | 2 and transmittances Ti = |si iOUt | 2 
(i = 2, 3) arc independent of r and invariant under the 
mapping (A,w — ojq) — > (—A,luq - co), which amounts to 
conjugating the expressions in curly brackets in Eqs. (14). 
Therefore it is sufficient to study the properties of R 
and Ti in the interval A e [0, ^n]. 

Apart from a global phase factor r, the device's spec- 
tral behaviour can be normalized to the mode coupling 
strength V and described entirely by three real param- 
eters: (1) the reduced operation frequency (lj — wq)/V, 
(2) the reduced cavity-port decay rate "f/V, and (3) the 
phase lag A + tt between the reflection r and trans- 
mission t of the direct waveguide-to- waveguide coupling. 
The role of the latter is better understood by remarking 
that in Eq. (14) cos A = |i|/(2|r|). This direct coupling 
term, when present, adds an extra interfering contribu- 
tion to both circulator transmittances, T 2 and T 3 . 

The spectral response of the cavity ports can be bet- 
ter understood by partial fraction decomposition of the 
second term between curly brackets in Eqs. (14). It is 
straightforward to show that 



S2,out = e 



2 cos A 



VI + 8 cos 2 A 3 

F(A) 



o 

_ gi arg(3 cos A+isin A) 



7e 



+in/3 



7 e 



0(A) 
-in/3 



7- i( 



-) 7-i(fc 



(15a) 



2 COS A 2 i arg (3 C os A+isin A) 

*3,out — e ^ r- . — — „ e 



VI + 8 cos 2 A 3 



7e" 



-in/3 



7e 



+ivr/3 



7 — i(u> — cv + ) 7 — i(uj — cj_) 



(15b) 



When the ports are perfectly decoupled (A = ±2 n )> tnc 
circulator transmission coefficients are sums of standard 
Lorcntzian lineshapes, L± 7r /3, centred at the (split) res- 
onant frequencies ui±. In all other cases the direct term 
-F(A) adds a continuum background coupling (with zero 
phase) and will thus in interaction with the Lorentzians 
add a sharp asymmetric Fano resonance to the circulator 
transmittance close to u>± . Fano resonances in integrated 
photonic structures have been predicted and observed in 
many configurations (see for instance Section V in the 
review paper Ref. 24). The additional presence of the 
coupling term V introduces an extra nonreciprocal asym- 
metry in the port transmittances. This results in the 
±i-7r phase shifts of the Lorentzians at the frequencies 
to± and in particular in the sign reversal of these phase 
shifts between the output ports. The impact of nonre- 
ciprocal coupling on Fano resonances in optical cavities 
is, up to our knowledge, studied here for the first time. 



C. Discussion: Ideal lossless case 

We shall now demonstrate the influence of these non- 
reciprocal Fano-type resonances on the isolation level of 
the circulator, defined as 



J = 10 log ^. 

J 3 



(16) 



Figure 2 shows the dependence of I on the reduced fre- 
quency shift (uj — ojq)/V and the reduced decay rate j/V 
for a few representative values of A. Several features can 
be distinguished in these plots. First of all, the spectral 
behaviour of the isolation level is strongly asymmetric 
around the resonant frequency ujq of the cavity, except 
for the limiting cases of zero (A = ^tt) and maximum 
(A = 0) direct coupling. This is obviously a signature 
of Fano interferences. Secondly, regardless of the level of 
direct coupling there always seems to exist a frequency 
for which infinite isolation (i.e. T 2 — or T3 = 0) can be 
obtained by properly adjusting the cavity-port coupling 
strength 7. In fact, it can be proven mathematically 
using Eqs. (14b) and (14c) that perfect isolation can in- 
deed occur at all values of A except ±|7r. Thirdly, the 
sign of the maximum (infinite) isolation level varies as a 
function of A. This implies that by "tuning" the level of 
direct waveguide-to- waveguide coupling the sense of the 
circulation can be reversed! Finally, the circulation also 
exhibits a spectral sign change as soon as direct coupling 
appears. As a result, the sense of circulation (clockwise 
or counterclockwise) changes outside a certain band of 
frequencies around the optimum frequency. This might 
have interesting applications for integrated optical add- 
drop functions. 

In practical applications one is also concerned about 
maximizing the bandwidth £?(/ m in), defined as the length 
of the frequency interval in which the magnitude of the 
isolation level exceeds a predetermined threshold I m i n . 
Figure 2 shows that this bandwidth is particularly large 
for vanishing (A = ^n) and maximum (A = 0) di- 
rect waveguide-to- waveguide coupling. In these cases one 
proves easily using Eqs. (14b) and (14c) that T 2 and T 3 
are even functions of (co — ujq) and, thus, have an ex- 
tremum at cj - Perfect isolation will therefore occur at 
lu = luq, provided that the cavity-port coupling level 7 
takes a specific value 700, equal to |V"|V3 for A = i-zr 
or to |y|/V3 for A = 0. Interestingly, however, the 
bandwidth B (/,„;„) for a given value of 7 min can some- 
times be increased by tuning 7 away from 7^. For in- 
stance, Fig. 2(d) shows that when A = 0, the 20-dB 
bandwidth is largest for 7 slightly smaller than 700. In 
fact, Eqs. (14b) and (14c) can be used to show that 
the maximum band widths £?(/ m in) obtainable in circu- 
lators with A = and A = ^tt are identical for all 
imin > 101og(31+8Vl5) w 18 dB. Nonetheless, the struc- 
ture with zero direct coupling (A = ^w) has the advan- 
tage of lesser sensitivity to perturbations of the cavity- 
port coupling level off its optimum value. In addition, for 
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FIG. 2 Isolation level I (in decibels) as a function of the reduced frequency shift (w — ujq)/V and the reduced decay rate 7/V 
for A = ^7r, ^7r, j-k and 0. The isolation contours are plotted every 5dB. 



this structure the value of 7 maximizing the bandwidth 
for any reasonably high value of I m i n is very close to 700 , 
the value at which perfect circulation becomes possible. 
This facilitates the device design. 

In Fig. 3 we have plotted the frequency dependence 
of the port transmittances separately for three scenarios: 
complete absence of direct coupling (A = ^ir), maxi- 
mum direct coupling (A = 0), and an intermediate level 
of direct coupling (A = jir). For each of the three 
scenarios, the port transmittances are plotted for two 
levels of the cavity-port coupling: 7 = 700 (top) and 
7 = j7oo (bottom). The features mentioned above ap- 
pear clearly. For instance, the sign change of the isolation 
level — T% > T3 T3 > T2 — when going from per- 
fectly decoupled to maximally coupled circulator ports is 
clearly visible. Likewise, as soon as direct coupling ap- 
pears (A 7^ ±^7r), one can observe how Fano resonances 
in the port transmittances are responsible for both the 
asymmetry and the spectral sign changes of the isolation 
levels. 



The appearance of the above features can be 
understood by considering the form of Eqs. (15) 
complex transmission coefficients. Neglecting the 
phase factor e lr , they contain two contributions: 
constant real term F(A) that varies from to | as 
from to and (2) the sum of two oppositely ^ir- 
shifted Lorentzians. The latter complex sum is 
plied by the factor —0(A) of constant magnitude, 



better 
for the 
global 
(1) a 
A goes 
-phase- 
multi- 
|, and 



phase that varies from ^ir to 0. As shown in Fig. 4, in 
the complex plane the Lorentzian parametric curves de- 
scribe circles with radius |, starting and ending at the 
origin and attaining their "apex" at e ±17r / 3 when u) = u>± , 
respectively. For the other port, the roles of the circles 
are interchanged. When these two parametric circular 
curves are added, one obtains the loci shown with the 
black curves in Fig. 4. In the absence of direct-pathway 
coupling, the Fano term (a) is zero and these curves, 
multiplied by the complex factor — 4>{\tt) — — |i, effec- 
tively describe the port transmission coefficients. For a 
suboptimal cavity-port coupling strength 7, at the cen- 
tral frequency oj — loq the values of the two Lorentzians 
contributing to S3 ;OU t [marked with triangles in Fig. 4(a)] 
do not exactly cancel each other. As a result, as shown 
in Fig. 5(a), complete isolation is not achieved. By ad- 
justing the value of 7 to 700 = Vy3, the speed with 
which the two circles are traversed can be tuned so that 
at lj = u>q the two Lorentzians take opposite, purely 
imaginary values, thus making S3 i0 ut vanish. This sit- 
uation is illustrated by Fig. 5(b). Fig. 5(c) shows a sim- 
ilar construction for an intermediate level of direct cou- 
pling (A = jn). This description makes the behaviour 
of the port transmittances for A ^ ±^ir easier to under- 
stand. The complex factor <fi(A) first "rotates" the com- 
plex paths L 7! / 3 (u — Lu±) + L_ 7r / 3 (uj — ujz i z) to the positions 
indicated by the dotted curves in Fig. 5(c). The Fano in- 
terference term F(A) subsequently shifts the paths along 
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-4 -2 2 4 -4 -2 2 4 -4 -2 2 4 

(lj — uo)/V (u — ujo)/V (u> — U}q)/V 



FIG. 3 Frequency dependence of the port transmittances following Eqs. (14b) and (14c) for (a)-(b) complete absence of 
direct pathway coupling (A = W); (c)-(d) an intermediate level of direct coupling (A = j7r); and (e)-(f) maximum direct 
pathway coupling (A = 0). For each considered A, the lower plot corresponds to optimal 7 = which gives either T2 or 
T3 = at a certain frequency. The upper one shows the behaviour when 7 is suboptimal, 7 = ^700- 




FIG. 4 Traces of the two complex Lorentzians occurring in (a) Eq. (15a) for S2, ou t and (b) Eq. (15b) for S3, ou t (with 
7 = -j7oo = t^) treated as complex functions of the frequency u. As oj goes from — 00 to 00, L n / 3 (w — oj+) and 
L-Tr/'ii 1 ^ ~ u -) follow the blue and red circular paths, respectively, in the direction of increasing line thickness, starting and 
ending at the origin. The black curve traces the sum of the two Lorentzians. The rate of change of the line thickness, as well 
as the locations of the markers corresponding to uuq and oj+, indicate the speed with which the curves are traversed. 



the real axis. The latter shift explains why the ratio 
log( I S2 1 / 1 S3 1 ) can change sign as a function of frequency: 
the paths do not start or end anymore in the origin of 
the complex plane. On the other hand, the asymmetric 
phase shift induced by 0(A) is responsible for the sharp 
asymmetries of the port transmittances. Only for max- 
imal direct coupling does the symmetry reappear, since 
arg(-0(O)) = 7T. 

In summary, an abstract coupled-mode-theoretical de- 
scription of a 3-port circulator including possible Fano 
resonances shows that the circulator ports do not neces- 



sarily have to be perfectly decoupled (as in Refs. 13, 14) in 
order to obtain reasonable behaviour. In fact, somewhat 
paradoxically, devices with maximum direct coupling also 
provide a large bandwidth, exhibiting, in addition, an in- 
teresting switching effect. This is a first indication that 
a technologically less demanding cavity layout not em- 
bedded in a PhC background — and hence without a 
guarantee of perfect port decoupling — is worth consid- 
ering. 
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FIG. 5 (a) Traces of S2, ou t and S3, ou t, treated as complex functions of the frequency u>, of a circulator with no direct-pathway 

port coupling (A = i7r), a suboptimal cavity-port coupling strength 7 = j7oc = ^-V and phase r = 0. Increasing line 
thickness indicates as before the evolution of uu from -co to 00. Note that because of the absence of direct coupling, the two 
paths are simply the black curves from Figs. 4(a) and (b) scaled by the complex factor 0(A = i 7r) = § e~" r/ ' 2 . Note also that 

Fig. 3(b) shows the squared modulus of these parametric complex curves, (b) Same for the optimum cavity-port coupling 
strength 7 = 700. The black curve is traversed twice. The squared moduli of these curves correspond to Fig. 3(a). (c) Same 
for A — j-k and 7 = 700 = ^ V. The dotted curves show Si, ou t after substracting the direct-coupling term F(0). Perfect 
circulation occurs now at a frequency different from uq. This case corresponds to Fig. 3(c). 



D. Inclusion of radiation loss 

Another, probably more detrimental, effect that can 
play an important part in the operation of circulators is 
the radiative energy loss. It can be fairly easily handled 
by the coupled-mode theory. We will limit our discus- 
sion to the case without direct coupling (A = ^7r), which 
in the previous subsection was found to yield maximum 
bandwidth and to be most tolerant to perturbations of 
the cavity-port coupling strength. We will also assume 
the energy to be lost solely owing to radiation of the cav- 
ity mode into free space, neglecting the energy radiated 
when the incident waveguide mode is reflected from the 
end of the input waveguide. This allows us to assume 
that the form of the coupling matrices C and D is unaf- 
fected. Note that Eq. (2a) remains valid if by T in this 
equation one implies only the decay rate to the waveguide 
ports. 

Following Ref. 25, pp. 208-212, we incorporate the cav- 
ity radiation loss in the coupled- wave equations (1) by 
replacing T by T + T r , where r r = 7r [ J ? ] and j r is 
the decay constant related to radiation loss. Proceeding 
analogously as before, we solve the coupled- wave equa- 
tions for the amplitudes of the outgoing modes in the 
three waveguides: 



sr.out — 



S2,out 



S3, out 



4 

+ r 



7[7 + 7r - K w 



3 [7 + 7r-i(w-w )] 2 + ^ 2 ' 

2 7^ + 7 r + V\/3 — i(w — wo)] 

3 [7 + 7r - \{u) - cu )} 2 + V 2 

2 7(7 + 7,. - VV3 - i(u - ujq)] 



3 [7 + 7r - i(w 



wo) 



V 2 



(17a) 
(17b) 
(17c) 



From eqs. (17b) and (17c) it can be immediately seen 
that despite the presence of loss, infinite isolation level / 
can still be obtained at frequency uj = luq provided that 
the coupling coefficient 7 is taken as 



7-1 



V\y/Z- 



7r- 



(18) 



This underlines that in the presence of radiation loss in 
the cavity the critical coupling of the waveguide ports 
to the cavity to obtain perfect isolation must become 
weaker. If the radiative decay rate 7 r becomes too high, 
the above condition cannot be fulfilled for any positive 
value of 7 and hence perfect circulation becomes impossi- 
ble. Similar conclusions can be reached if the waveguide 
ports are not perfectly decoupled (A 7^ ^ir). Interest- 
ingly, for the case of A = ^tt, when 7 is chosen along 
(18), the isolation level I(uj) [and the bandwidth B(I m i n )] 
becomes independent from 7 r : 



I{uj) = 10 log 



1 



12V 2 



(U - LOof 



and 



-B(^min) 



4\/3|V| 



VlQWio _ 1 ■ 



(19a) 



(19b) 



However, the power transmitted to the 'hot' waveguide, 
T 3 , obviously decreases when the loss coefficient j r in- 
creases. With 7 given by Eq. (18), the expression for the 
maximum value of T3 takes the form 



T 3 (u = uj ) 



(V3-Jr/\V\) 2 



(20) 



This formula is important because it allows to deter- 
mine the minimum quality factor Q T = o;o/(27 r ) of a 
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cavity with a given relative frequency splitting Aw/wo 
that can be used to build a circulator with a prescribed 
minimum level of transmitted power, T3 im ; n . (Note that 
Q T describes solely the decay caused by radiation to free 
space, rather than that due to interaction with waveguide 
modes in the circulator ports.) Indeed, expressing j r in 
Eq. (20) in terms of Q r and using the relation Aw = 2\V\, 
we obtain that T 3 (w = uj ) > T min if and only if 



Qr > 



— (1-VT~)V3 



(21) 



Thus, a circulator built with a cavity having Aw/wo = 
0.001 (a typical value) can provide 50-percent peak trans- 
mission if Q r > 1970. A quality factor Q r > 11,250 is 
needed for 90-percent efficiency, and Q r > 115,000 for 
99-percent. 



III. RIB- WAVEGUIDE-BASED CIRCULATORS 

A. Design of axisymmetric ring cavities 

As is clear by now, the performance of a cavity-based 
circulator depends heavily on the strength V of the MO 
coupling of the even and odd cavity modes. As can be 
seen from Eqs. (18) and (19), an increase of \V\ augments 
the bandwidth, necessitates less critical coupling and 
leads to higher tolerance towards radiation loss. Wang 
and Fan [23] derived a pcrturbational formula for the 
mode coupling strength V of a cavity containing a MO 
material polarized in the z direction. Such a material is 
characterized by the permittivity tensor 



e ig 
-ig e 
e 







(22) 



where g, proportional to the material's magnetization, is 
the gyroelectric index. From now on we shall assume g 
to be real. Together with the real-valuedness of e this is a 
necessary condition for a Hcrmitian e and thus a lossless 
MO material. Wang and Fan [23] obtained [30] 



V = iV = - 



i cjq J g(r) z ■ [E* e (r) x E {r)] dr 
2 ^Je(r)\E (r)\ 2 drJe(r)\E c (r)\ 2 dr 



where the integrals run over the whole cavity volume, 
E c and E are the electric fields of the even and odd 
modes of the non-magnetized cavity, and uj is their fre- 
quency Thus, V is proportional to the cross product of 
E* and E weighted by the off-diagonal component of 
the permittivity tensor, ig. 

It is convenient to introduce the dimensionless reduced 
MO coupling strength v defined as |V"|/(wo|<7|max), where 
Mmax denotes the maximum magnitude of g over the 
cavity volume. The parameter v depends then solely on 
the geometry of the cavity and not on the strength of 



the MO effect in the chosen material. It can therefore 
be used to compare the merits of cavities with different 
geometries. 

By exploiting the analytical properties of the eigen- 
modes of the Helmholtz equation in axisymmetric ge- 
ometries, we will now derive a general design strategy 
for resonant 2D cavities with large frequency splitting 
in a uniform static external magnetic field. In contrast 
to general PC cavities, circularly symmetric cavities can 
be handled analytically, which makes it possible to get a 
better insight into their properties. 

We consider, then, a system of concentric rings [see 
Fig. 6(a)], described by a piecewise-constant relative per- 
mittivity e(p) independent from the azimuthal coordi- 
nate <j) an d the vertical coordinate z. The relative per- 
meability p is taken to be 1 everywhere. In the p- 
polarization case [waves propagating in the (p, <fr) plane 
with the magnetic field parallel to the z axis] Maxwell's 
equations 



ldH z 



dH z 
dp 



—iujee E p , (24a) 
iojeeoE^, (24b) 



d -^ + -E,-- d -^=^H z (24c) 
op p p dip 

reduce to the Helmholtz equation for the z component of 
the magnetic field, H z : 



d ndH z \ l dH z 1 d 2 H z tj*_ _ Q 
dp\e dp J p dp p 2 dif> 2 c 2 z 



Here ui denotes the frequency, eo and po the permittiv- 
ity and permeability of free space, and c = 1/ ^e p the 
speed of light in vacuum. By the usual method of sepa- 
ration of variables one can show that the eigenmodes of 
the system occur in degenerate pairs of "even" and "odd" 
modes with magnetic field of the general form 

H ze (p, <t>) = Ri(p) cos(/0) and H zo (p, <j>) = R t (p) sin(^), 

(26) 

respectively, where the azimuthal order I is an integer 
and the radial dependence Ri(p) is a solution of 



e(p) 



d_ 

dp 



1 dR l 
e(p) dp 



I dRi 
P dp 



+e(p)k% 



1- 



l 2 



e(p)k 2 p 2 



Ri(p) = o, 

(27) 



which reduces to the Bcsscl equation within each uniform 
ring. In the above equation ko = lu/c is the free-space 
wave number. Within each ring, Ri(p) can then be ex- 
pressed as a superposition of the Bessel functions of the 
first and second kind: 

Ri(p) = a m Ji{n m k(ip) + b m Yi(n m k p), (28) 

where m is the ring's number, n. m (= y r e^) its refractive 
index, and a m and b m constant coefficients. 

From Eq. (26) and Maxwell's equations (24) we can 
derive the general expressions for the cross product of 
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0.0 0.5 1.0 1.5 2.0 2.5 

P (pm) 

FIG. 6 (a) Magnetic field H z of the even eigenmode of 
azimuthal order I — 1 and wavelength A = 1.300 pm 
supported by the cavity composed of an alternating sequence 
of rings made of air and a MO material with refractive index 
2.25. (b) Radial dependence of the magnetic field H z of the 
cavity eigenmodes and of the cross product of their electric 
fields. Both functions are normalized to their maximum 
values. The areas filled with the MO material are shaded. 



the electric fields corresponding to a pair of degenerate 
modes with frequency ujq and azimuthal order I, 



z-(E* c x E Q ) 



and for their squared norm, 



lRi(p)Bl(p) 

"Ip4 Hp)? 1 



(29) 



J t{r)\E c {r)\ 2 dr = J e{r)\E {r)\ 2 dr 



u 2 p4 6 (p) 



■ dp, 
(30) 



Clearly, the sign of the integrand in the numerator of 
the above expression depends on the sign of the product 
g(p)Ri(p)R' l (p). In general, this sign will oscillate due 
to the oscillatory behaviour of the function Ri(p) itself. 
Therefore, to avoid unnecessary cancellations in the in- 
tegral in question, the ring boundaries need to be placed 
so that the product g(p)Ri(p)R[(p) be always nonneg- 
ative (or nonpositive). Short of introducing inversely- 
magnetized domains, this can only be achieved by ensur- 
ing that the boundaries between the MO and non-MO ma- 
terial coincide with the zeros of the product Ri(p)R[(p), 
i.e., the zeros and extrema of Ri(p). With this design 
principle in mind, we offer the following method for the 
determination of the radii of the alternating MO and non- 
MO rings making up a cavity supporting a pair of degen- 
erate modes with azimuthal order I > 1 at a particular 
frequency ujq = k^c. 

As a first step, we determine the radius po of the cen- 
tral rod ("zeroth ring"), assumed to have refractive in- 
dex uq. In this rod, the radial field dependence is ex- 
pressed solely by means of the Bessel function of the first 
kind Ji(nokop), since Yiin^kop) has a singularity at the 
origin. Hence, bo = and, since the global mode am- 
plitude is arbitrary, oq can be set to 1. It is then evi- 
dent that the product Ri(p)R[(p) in the central rod will 
not change sign until the first extremum of the Ji(nok p) 
function. Therefore we set po to the value of p at which 
this extremum occurs. 

The amplitudes a\ and b\ of the Bessel functions in 
the first ring, having refractive index n\, can now be de- 
termined from the condition of continuity of H z and E$. 
The outer radius of the first ring, p\ , should then be cho- 
sen so as to coincide with the first zero of the function 
aiJi(riikQp) + 6 1 Y/(n 1 /c p) located in that ring: this is 
where the product Ri(p)R' l (p) will again change sign. By 
repeating this process, we obtain the radii of the subse- 
quent rings. The outer radii of the even rings (with re- 
fractive index no) coincide with the extrema of the func- 
tion Ri(p), and those of the odd rings (with refractive 
index rii), with its zeros. Thus, the cavity designed in 
this way resembles an annular Bragg grating [21], with 
each layer a quarter- wavelength thick (in the sense of the 
Bessel- function quasi-periodicity) . 

We need now to make sure that the field of the mode 
constructed in this way is localized, i.e., that Ri (p) decays 
as p tends to infinity. Let us estimate the ratio of Ri (p) 
at its two successive extrema, located at p2n and /?2n+2, 
where n 3> 1. The Bessel functions of large arguments 
can be approximated by [26, Eqs. (9.2.1) and (9.2.2)] 



with R[{p) denoting the derivative of Ri(p). Substituting 
Eqs. (29) and (30) to the general formula for the coupling 
strength (23), we obtain 



Jl{x) 
Yi(x) 



2 I ( <»s 

irx I sin 



7T 7T 

x I 

2 4 



(32) 



V 



. , , g(p)Ri(p)R[(p) J 
1W ° ' 7{pf p 



PIMp)} 2 + p 2 [R[{ P )] 
Pe(p) 



dp 



(31) 



By construction, the field in the (2n+l)th ring, located 
between pm and p2n+i, behaves approximately as 



Ri(p) ~ c 2n +iP 1/2 cos[ni& (p - P2n)], 



(33) 
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and in the (2n + 2)th ring as 

Rl(p) ~ c 2n+2 p~ 1/2 sin[n fc (p - p 2n +i)}, (34) 

where c 2n +i and c 2n+2 are constant coefficients. Since the 
rings are supposed to be a quarter-wavelength thick, the 
expressions nik (p 2n+ i - p 2n ) and n k (p 2n+2 - p 2n +i) 
are approximately equal to ^ and the successive extremal 
values of Ri (p) are 



Rl(p2n) « C 2n+ i p 2 n 2 , 



(35a) 



Rl(p2n+2) ~ C 2n+ 2 P2n+2 ~ c 2n+2 









P2 " + i 


\n 


>-)] 



(35b) 



From the requirement of continuity of E$ oc R[(p) /e(p) at 
P = P2n+i it follows that c 2n+ 2/ c 2n+i = -n /«i- Hence, 



Rl(p2n+2) 



Rl{p2n) 



no 
ni 



2k P2r, 



1 1 
1 

n ni 



'1/2 



n 



the cavity. For systems based on garnets, this approach 
might not be practical, since it would be necessary to 
fill the (thin!) slits between MO rings with a material 
of refractive index even higher than ~2.5. On the other 
hand, this approach would be very well suited to systems 
based on low-index MO polymers, such as those analysed 
recently by Jalas et al. [6]. 



B. Design of ring-cavity circulators: introduction 



-1/2 



m 
(36) 

Thus, a localized mode is obtained only if the refractive 
index of the central rod, no, is chosen smaller than that 
of the first ring, n\. 

Figure 6(a) shows the geometry of an example cavity 
constructed with the above algorithm. The system con- 
sists of an alternating sequence of rings made of air (re- 
fractive index no = 1) and a MO material with refractive 
index n\ = 2.25. This value was chosen as the effec- 
tive index of the fundamental s-polarized guided mode 
of the heterostructure composed of an air-covered 340- 
nm-thick layer of BIG (refractive index 2.51) grown on a 
galium-gadolinium-garnct (GGG, refractive index 1.97) 
substrate, at wavelength A = 1.3 pm. The radii of the 
first seven rings are listed in Table I. This cavity supports 
a pair of degenerate localized eigenmodes of azimuthal or- 
der 1 = 1. The magnetic field of the even mode is plotted 
in Fig. 6(a); the field of the odd mode can be obtained 
by rotating the map from Fig. 6(a) by 90° in the counter- 
clockwise direction. In Fig. 6(b), the radial dependence 
of these magnetic fields, Ri(p), is juxtaposed with that of 
z-El x E , calculated from Eq. (29). Clearly, Ri{p)R[{p) 
is negative throughout the MO rings and positive else- 
where; as a result, no cancellations in the upper integral 
in Eq. (31) occur. The reduced coupling strength of this 
pair of modes is v — 0.00874. This is more than an or- 
der of magnitude stronger than the corresponding value 
for the cavity proposed in Ref. 13 if it were placed in a 
uniform magnetic field: v = 0.0006. 

The plot in Fig. 6(b) shows that z-E* c ~xE a takes much 
larger values in the even-numbered (low-index) rings than 
in the odd-numbered (high-index) ones. This is because 
this cross product, given by Eq. (29), is inversely propor- 
tional toe 2 . A potential way to increase the reduced cou- 
pling strength v still further consists therefore in placing 
the MO material in the even- numbered rings. Of course, 
owing to the constraint (36), the MO material would in 
this case need to become the lower-index component of 



The next step in the circulator design consists in cou- 
pling the MO ring cavity with waveguides. Side coupling, 
which is often used for classical ring resonators (consist- 
ing of only one ring) [27], is not applicable in our case. 
As Fig. 6(a) shows, far from the centre of the cavity, 
the field of the cavity mode varies very slowly in the az- 
imuthal direction, so it could only be phase-matched by 
a waveguide mode with an effective index much lower 
than the refractive index of air — which cannot be ob- 
tained in practice. It is therefore preferable to transfer 
energy between the cavity and the waveguides by means 
of butt-coupling. 

Figure 1 shows the general geometry of the class of 
structures we have considered. They consist of a resonant 
cavity, composed of ni full and n s split rings with inner 
and outer radii determined by the procedure described 
in the previous subsection, and three radially oriented 
identical waveguides of width d wg . The distance from 
the centre of the cavity to the ends of the waveguides is 
p wg , while the widths of the slits in the split rings are 
denoted by d<f, n with n = rif + 1, rif + 2, . . . , nf + n s . The 
rings and waveguides are made of a MO material with 
permittivity 



£h iff 
-iff £h 
e h 



(37) 



and are embedded in an isotropic medium with permit- 
tivity e\. 

In our simulations, the permittivities £h and ei were 
chosen as in the previous subsection, i.e. as £h = (2.25) 2 
and ei = l. A Faraday rotation of 10-20°/pm has been 
measured in BIG at A = 620 nm [19]. The tail of the 
Faraday spectrum in this reference approaches a couple 
of deg/pm, which corresponds to g ~ 0.05-0.1, close to 
the near-infrared part of the spectrum. To assess the 
maximum performance that can in principle be achieved, 
we fixed g = 0.1. 

The level of coupling between the cavity and waveguide 
modes, and hence the circulator's performance, depends 
of course on the values of all the geometrical parameters, 
which should therefore be optimized. The space spanned 
by them is rather large, and it is not possible to scan it 
exhaustively. Therefore our optimization of the presented 
structure has been somewhat heuristic. The waveguide 
width d wg was fixed to 250 nm. The radii of the MO 



Ring number Inner radius (nm) Outer radius (nm) 



1 


381 


539 


2 


847 


998 


3 


1309 


1457 


4 


1772 


1919 


5 


2236 


2382 


6 


2700 


2846 


7 


3165 


3310 



TABLE I Radii of the high-index rings of the cavity with 
eh = (2.25) 2 and ei = 1 supporting a resonant mode with 
azimuthal order I = 1 at wavelength A = 1300 nm. 



rings, listed in table I, were determined with the proce- 
dure described in the previous subsection to ensure the 
existence of a pair of cavity modes with the azimuthal 
order I = 1 at the wavelength A = 1300 nm. 

With the chosen value of g, the relative frequency split- 
ting of these modes is Aw/wo = 2gv = 0.00175. From 
Eq. (18), in the absence of losses the optimum value of 
j/loo is 0.00151, i.e., the quality factor Q = cj /(2j) de- 
scribing the cavity-waveguide coupling should be Q = 
330. The quality factors of the chosen cavity with 3 and 
4 rings are 163 and 829, respectively [31]. Therefore one 
can expect that the waveguides of an optimally designed 
circulator should end somewhere close to the third inner- 
most ring — or possibly even inside it, since the coupling 
to waveguides is doubtlessly less efficient than that to the 
whole surrounding free space. 

The total number of rings necessary for ensuring a pre- 
scribed level of peak transmission T m ; n could in principle 
be estimated from Eq. (21): for instance, for T m ; n = 0.9 
the quality factor Q r describing radiation loss should be 
greater than 6430. This is already ensured by a 6-ring 
cavity, whose quality factor reaches 21,140. However, the 
quality factors of cavities with split outer rings will nec- 
essarily be smaller than of those with full rings; therefore, 
a larger number of rings might be necessary to ensure a 
90-percent efficiency. In our calculations, we restricted 
our attention to systems with at most 7 rings. 

Having fixed the number of split and full rings of the 
circulator in this way, we were essentially left with the 
problem of optimizing the values of d$ n and p wg . The 
ultimate figure of merit, the bandwidth, is nonzero only 
for structures already rather close to optimum, which 
makes it a cumbersome objective function. Therefore we 
optimized instead the maximum isolation level, hoping 
that circulators with large values of this parameter would 
also be characterized by a large bandwidth. 

Numerical calculations of the transmission through the 
circulators studied in this section were performed with 
the finite-element method (FEM) using the RF module 
of the COMSOL software package. To account for the 
presence of idealized semi-infinite waveguides, the com- 
putational domain, shown in Fig. 7, was constructed as 
follows. The region surrounded by the dashed line is 
a fragment of the physical system shown in Fig. 1. The 
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FIG. 7 Geometry of the domain used in FEM calculations. 

parts of its boundary lying "far" from the waveguide ends 
are covered with perfectly-matched layers (PMLs) imple- 
mented by means of the complex coordinate transform 
P !-> spML(l+i)(p— Pin)/rfpML, where p is the radial coor- 
dinate measured form the centre of the cavity, p m denotes 
the radius of the inner PML boundary, c?pml = 600 nm 
is the PML thickness, and spml = A its strength. The 
radius p ln is chosen so that the distance from the outer- 
most ring to the inner PML boundary is d scp — 5000 nm. 
Perfect-electric-conductor boundary conditions are im- 
posed on the parts of the domain's boundary adjacent 
to the PML. In constrast, the electromagnetic fields on 
the segments P n (n = 1, 2, 3) of length d port = 2250 nm, 
perpendicular to the waveguides, are constrained to be a 
superposition of the incoming and outgoing guided modes 
of the corresponding waveguides. The profile of these 
modes is calculated analytically and normalized to uni- 
tary power, and the amplitude of the incoming mode is 
set to unity on Pi and to zero on P 2 and P 3 . Phys- 
ically, these constraints correspond to the assumption 
that all the radiative waveguide modes excited by the 
cavity decay before reaching the ports P n . The domain 
is divided into triangular nodal Lagrangian elements of 
order p = 5 and maximum allowed size /imax/v^ where 
ft max = 500 nm and e is the local permittivity. 

We have verified that above choice of the com- 
putational parameters allows to calculate the 20-dB- 
bandwidth of typical circulators with an accuracy of 
about 3% (see Rcf. 28, Section 4.5.2, for a detailed con- 
vergence study). 

C. Design of ring-cavity circulators: geometry 
optimization 

We initially focused on circulators with d^ n indepen- 
dent from n, i.e., with waveguides enclosed in rectangular 
"slits" of width d^. The contours in Figure 8(a) show the 
geometry of the best structure we have found, called CI 
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in the following. It has 3 full and 4 split rings. The slit 
width g?0 = 1770 nm and the waveguide ends are located 
at p wg = 1240 nm from the cavity centre, so that the 
waveguides cross the outermost full ring. For ease of ref- 
erence, the values of all the geometrical parameters of CI 
are listed in table II. This circulator offers maximum iso- 
lation level of 35 dB, and the wavelength dependences of 
T 2 , T 3 and I are shown in Fig. 8(b). Clearly, the curves 
are fairly symmetric with respect to the central wave- 
length A = 1300.0 nm, which indicates that the direct 
coupling between waveguides is insignificant. Figure 8(c) 
shows the map of the magnetic field at A = 1300.0 nm. 
At this wavelength, 88% of the input power is transmit- 
ted to waveguide 2; the rest is mainly lost to the sur- 
rounding free space. Far from the peak, the amount of 
these losses can exceed 50%. This behaviour contrasts 
with that of the PC circulator, where almost 100% of 
the input power remains in the waveguide system due to 
the quasi-perfect isolation provided by the surrounding 
periodic lattice. The 20-dB bandwidth _B(100) of circu- 
lator CI is 0.729 nm (129 GHz). 

The grey curves in Fig. 8(b) show the transmittance 
curves of this circulator as predicted by the coupled- wave 
model with radiation losses taken into account, presented 
in Section II D. The values of the parameters wo, Acj, 7 
and 7 r were found by fitting the expressions for T2 and T3 
obtained from eqs. (17) to the values calculated numer- 
ically. The Levenberg-Marquardt algorithm was used as 
the fitting procedure. The best fit was obtained for pa- 
rameters corresponding to Ao = 2ttc/ljq — 1299.9 nm, 
AA = 27rcAw/cjg = 2.3 nm, Q = u /{2l) = 370 and 
Q r = ujo/(2-f r ) = 5730. Clearly, there is a good match 
between the theoretical and numerical curves; its qual- 
ity would probably be further improved by taking into 
account the direct coupling between waveguides, which 
causes the slight asymmetry of the numerical plots. The 
quality factor related to losses, Q r , is significantly lower 
than that of an isolated cavity with 7 full rings, which is 
as large as 107,000. This is obviously due to the presence 
of slits. On the other hand, the position of the ends of 
the waveguides (just inside the third innermost ring) is in 
good accord with the predictions made in Section IIIB. 

We have found this device very tolerant to variations of 
the slit width d<p; Figs. 9(a)-(b) show the dependence of 
the maximum isolation level and bandwidth of CI on this 
parameter. It can be seen that the bandwidth stays above 
0.5 nm in a 300-nm-wide range of d<f>. The constraints 
on p wg are more stringent: as shown in Fig. 9(c), the 
corresponding range of p wg is about 40-nm wide. 
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TABLE II Geometrical parameters of the circulators 
analysed in the text. The ring radii are listed in the first 
(rif + n s ) rows of table I. 
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FIG. 8 (a) Magnitude of the magnetic field in circulator CI 
at the wavelength A = 1300.0 nm, corresponding to the 
maximum isolation level. The waveguide mode is incident 
from the left, (b) Wavelength dependence of the 
transmission (left axis) and isolation level (right axis) of CI. 



The performance of this structure changes rather 
abruptly when the number of rings is modified. For in- 
stance, if the seventh ring is removed, the maximum iso- 
lation level decreases to only 21 dB and the bandwidth 
to 0.386 nm (69 GHz). However, it is possible to improve 
these figures by readjusting the slit widths and the po- 
sition of the waveguide ends: for p wg — 1250 nm and 
c?0 = 1570 nm (10% less than in the 7-ring case) I max 
reaches 33 dB and the bandwidth 0.603 nm (107 GHz). 

Having noted that both for 6 and 7 rings the optimum 
angular length of the removed sectors of the outermost 
ring is almost the same, ^34°, we test the performance 
of a second class of structures, in which the outer rings 
are truncated along radial lines instead of ones parallel 
to the waveguides. The latter are thus enclosed by coni- 
cal rather than rectangular air slits, as illustrated in the 
inset of Fig. 10. We found the optimum cone aperture 
0conc to be 35.1°, close to the value cited above. The 
optimum position of the waveguide end, p wg = 1210 nm, 
is also only slightly different from the original one. More 
importantly, the maximum isolation level and bandwidth 
decrease much less (from 35 to 28 dB and from 0.770 to 
0.708 nm, or 137 to 126 GHz, respectively) when the sev- 
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FIG. 9 Tolerance of the bandwidth B(20dB) and the 
maximum isolation level J ma x of circulator CI to 
perturbations of the parameters (a)-(b) and (c)-(d) p wg . 
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FIG. 10 Wavelength dependence of the transmission (left 
axis) and isolation level (right axis) of circulator C2 (inset) 
with waveguides enclosed in conical slits. 



enth ring is removed. This relative insensitivity to the 
details of the geometrical structure of the exterior region 
of the device is the behaviour that one would intuitively 
expect from a well-designed circulator; therefore, struc- 
tures with conical slits seem closer to the ideal than those 
with rectangular ones. Figure 10 shows the wavelength 
dependence of T 2 , T 3 and I of the optimum 7-ring circu- 
lator with conical slits, referred to as C2. Its parameters 
arc listed in table II. 



IV. PERSPECTIVES — TOWARDS 3D 
NONRECIPROCAL CAVITIES 

A. High-contrast 3D cavities 

It is well known that the properties of a 3D cavity 
depend crucially on its vertical multilayer profile, and a 
cavity with a high Q in 2D can leak heavily in 3D if the 
presence of the substrate and superstrate is not taken 
into account. The design algorithm of Section III A is 



based on 2D considerations, and a similar algorithm in 
3D is nearly impossible to develop. However, the 2D al- 
gorithm has allowed us to compare our designs to previ- 
ously reported circulator cavities [13, 23, 29]. So far, the 
true 3D nature of a real integrated system has only been 
taken into account in our work by using an effective-index 
approximation for the unetched MO rings. The material 
between the rings has, however, always been set to be air 
— an approach that is reasonable only if the garnet layer 
is entirely etched away between the MO rings. Such a 
high-contrast ring cavity cannot, however, be fabricated 
in a membrane form, since the rings are not intercon- 
nected. Therefore, in practice the rings will need to be 
placed on a solid substrate with refractive index rt su b- 
Far from its centre, the ring cavity can be locally ap- 
proximated by a linear grating with period 



A = - — + — 

4 \n nx 



(38) 



If A > A/2n su b, all Bloch eigenmodes of the grating will 
leak into the substrate. In order to avoid that, n su b must 
be smaller than A/2A. For no = 1, m = 2.25 and A = 
1300 nm, the substrate index must be less than 1.4. This 
might prove very difficult to achieve. 



B. Low-contrast 3D cavities 

In low-index-contrast cavities substrate leakage is 
much easier to avoid, albeit at the price of larger cavity 
due to smaller Bragg gap. This can be seen by rewriting 
the no-leakage condition as n su b/2 < n e ff j o^eff,i/("efF,o + 
n e s,i)- This condition is always met if both effective in- 
dices are larger than the substrate index n sub , around 1.9 
for typical garnet substrates (such as GGG). 

We have considered two classes of low-contrast cavi- 
ties. In one of them, shown in Fig. 11(a) and labelled 
Gl henceforth, the garnet layer (cfeiG = 280 nm) is left 
unetched and the cavity is formed by etching in a high- 
index cladding layer (n sp i = 3.5, d sp \ = 80 nm). The 
MO properties of the exposed BIG regions are assumed 
to be destroyed by e.g. ion implantation. In the second 
class of cavities, shown in Fig. 11(b) and referred to as 
G2, rings etched in BIG (c?big = 330 nm) are buried in a 
material with index n sup approximately equal to that of 
the GGG substrate. A potential candidate is silicon ni- 
tride (SisN4); around A = 1300 nm, the refractive index 
of both GGG and S13N4 is 1.97. As indicated above, the 
use of the effective- index approximation to determine the 
position and width of the rings by help of the algorithm 
of Section III A is now better justifiable. We have used 
an in-house-developed rigorous FEM code [28, Section 
5.3] to compute the resonant frequencies of axisymmet- 
ric 3D cavities belonging to the just described classes, 
crosschecking the data obtained with the 2D model. Ta- 
ble III summarizes the results of these calculations. 

They prove first of all that the 2D design procedure 
is flexible enough to let us adapt it to 3D layouts. The 
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FIG. 11 Schematics of 3D structures of types (a) Gl and 
(b) G2. 
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a After numerical optimization of ring radii. 

TABLE III Relative frequency splitting Aw/wo and quality 
factor Q of the eigenmodes of several Gl- and G2-type 
20-ring cavities designed to support modes of azimuthal 
order ±L The values obtained with rigorous 3D FEM 
calculations are compared to the estimates obtained with 2D 

calculations in the framework of the effective-index 
approximation. Note that Q is defined as the average of the 
quality factors of the modes with aximuthal order I and — I. 



frequency splittings calculated with the FEM code are 
only slightly reduced with respect to those predicted by 
the 2D model. In fact this decrease can be entirely at- 
tributed to a reduced mode confinement in the MO layer 
due to the vertical permittivity profile. The more wor- 
rying aspect, however, is the strong diminution of the 
quality factor. The quality factors of the 3D cavities are 
an order of magnitude smaller than those obtained with 
2D calculations (for the same number of rings and related 
only to the in-plane confinement). 

The negative impact of large radiation losses (and thus 
low Q factors) on the operation of the nonreciprocal cav- 
ity can be easily understood by considering the following 
equations of the coupled-mode model including radiation 
losses: 



S3,out 

1 1 



2 7 [ 7 + 7r - VV3 - i(u - uj ) 
3 [7 + 7r -i(u;-u;o)] 2 + y 2 
Alo 



w 



(39a) 



(39b) 



where Q and Q r describe respectively the "useful" cavity- 
waveguide coupling and the radiation loss. The latter 
equation is obtained by rewriting Eq. (18) using 2\V\ = 
Alj and the definition of the quality factor Q = For 

simplicity only the case without direct coupling (A = ^ir) 
is considered. Equation (39b) expresses the optimum re- 



lationship between the relative frequency splitting and 
the two Q factors. It implies that for a radiative quality 
factor Q T lower than ujo / '(Aoj\/3) perfect circulation can- 
not occur, since it would require a negative (unphysical) 
value of Q. For such small Q r the coupling necessarily be- 
comes overcritical, i.e. j r + 7 is larger than the optimum 
value. 

Fig. 2(a) shows that the circulator has some toler- 
ance with respect to 7. At a ten-percent overtuning 
(7 <~ 2\V\) maximum isolation level of ~20 dB is still pos- 
sible. This will of course come at the price of a heavily 
reduced transmission level. Indeed, assuming a modest 
overtuning (in order to keep the isolation level reason- 
able), from Eq. (39a) we obtain maxT 3 = max|s 3;OUt | 2 ~ 
(-^) 2 , with 7 = |q describing the cavity- waveguide 
coupling alone. Thus, a ten-percent overtuning with 
Q r ~ uj /(AujV3) and Q ~ 0.1w /(Aw\/3) would lead 
to transmission levels of about 1%. 



C. Towards low-contrast nonreciprocal cavities 
with higher quality factor 

Table III shows that none of the 3D cavities under 
consideration supports a dipolar (1 = 1) mode with an 
inverse radiation quality factor 1/Q r smaller than 400% 
of the optimum coupling strength, Auj^/3/uj - The most 
direct way to increase Q r consists in moving towards de- 
signs for higher azimuthal orders I. The stronger az- 
imuthal variation leads to higher in-plane ^-components 
in the cavity mode and thus lower substrate leakage. 
Conveniently, according to (31), the relative frequency 
splitting also increases with azimuthal order. The price 
to pay is a higher modal volume. Table III shows never- 
theless that the increase in Q r is limited. The best cavity 
designed for I = 10, the one belonging to the class G2, 
supports a mode with 1/Q r of about 70% of the optimum 
coupling strength Auj^/Z/luq. As a result, at optimum 
coupling, 7 = f± ~ ^ 0.3^^ = 0.3\/3|y| and one 
still cannot expect transmission levels larger than 10% in 
the best-case scenario. 

We have recently applied a numerical optimization pro- 
cedure to adjust the positions and widths of the rings. In 
a first attempt this has already led to an improved de- 
sign for a 20-ring G2-type cavity (for an I = 10 mode). 
Fig. 12 shows the geometry of the original and the op- 
timized G2-type cavity, as well as the maps of the nor- 
malized magnetic field of their modes. The Q factor of 
the mode of the optimized cavity is as large as 2420: 
an increase by almost an order of magnitude from 334. 
Moreover, it is higher than that of the mode of the orig- 
inal 2D 20-ring cavity, 2139. This means that the Q 
factor of the optimized cavity is constrained by in-plane 
rather than out-of-plane radiation loss, and so it could be 
further augmented if necessary by increasing the num- 
ber of rings. At the same time, the relative frequency 
splitting stays quite tolerant with respect to this opti- 



16 



1 -| 

S 

3 A 



-1 - 



(a) \H Z \, original cavity 
CDDDDDDDDDDDDDDD 



3 



(b) \H Z \, optimized cavity 

EDDDnnnDEmDEmnEiDn 



| 
i 



1.0 



0.5 



0.0 
1.0 



0.5 



0.0 

2 4 6 8 

p (H 

FIG. 12 Magnitude of the z component of the magnetic 
field of the eigenmode of (a) the original G2-type cavity 
designed for I = 10, and of (b) its numerically optimized 
version; both in the absence of magnetization. The 
rectangles denote the cross-sections of BIG rings immersed 
in a material with refractive index 1.97. The improved field 
confinement (and thus higher Q) is easily seen. The MO 
interaction described by the cross product E* x E a (not 
shown here) is almost not influenced. 



mization, indicating that as far as MO effect goes the 2D 
design algorithm proposed here is already quasi-optimal. 
The obtained Aw/wo is 0.00261, even slightly larger than 
in the original structure. Using similar estimations as 
above, these values imply that in a best-case scenario 
perfect circulation with a maximal transmission close to 
90% is possible. The total device footprint remains close 
to (10A) 2 . This work together with a 3D analysis of the 
waveguide-cavity coupling problem will be the subject of 
a future communication. 



V. CONCLUSION 

In this work we have presented an extensive study of 
a new concept for a very compact and technologically 
feasible type of integrated optical circulator. Its opera- 
tion is based on a MO frequency splitting of a pair of 
degenerate cavity modes of opposite (azimuthal) parity. 
In contrast to previous reports on MO cavity circulators, 
the type proposed and studied here operates without op- 
positely magnetized domains or a full fc-space bandgap. 
The confinement is provided by a radial Bragg mirror 
consisting of uniformly magnetized, concentric MO rings. 



An original design algorithm taking into account the an- 
alytical properties of the functions describing the radial 
profile of the cavity modes allows to determine the opti- 
mum position and width of the rings. We have extended 
an abstract coupled-mode model to verify the impact of 
radiation losses and direct transmissions between the cir- 
culator ports. These might arise due to the absence of a 
full bandgap isolating the cavity ports. 

The extended coupled-mode model has indicated that 
the presence of direct pathway transmissions leads unex- 
pectedly to interesting alternative operation regimes of 
these nonreciprocal cavities. The combination of non- 
reciprocity and Fano-type interactions allows to achieve 
circulation regimes where the sign of the circulation can 
be spectrally tuned. Howsoever interesting the regime at 
high direct coupling is, devices designed to operate with- 
out direct coupling are expected to be less sensitive to 
perturbations of the cavity-port coupling level from its 
optimum value. We have demonstrated numerically the 
satisfactory circulator operation of a 6- and 7-ring cavity 
partially intersected by three 120°-spaced rib waveguides. 
Circulation extinction as high as 35 dB, and an inser- 
tion loss below 1 dB over a 20 dB circulation bandwidth 
of 130 GHz (around 1300 nm), have been demonstrated. 
The device is found to have tolerances on the position 
and dimensions of its features that are well within the 
limits of modern nanotcchnology. 

Looking forward to realistic 3D device layouts, we have 
studied different types of high- and low-contrast cavities. 
An FEM tool for the calculation of the eigenmodes of 
axisymmetric, 3D MO cavities has been developed for 
that purpose. Air-suspended high-contrast ring cavities 
being physically impossible, the minimum requirements 
for the Q-factor of the unavoidable low-contrast cavities 
have been indicated. First attempts at geometry opti- 
mization of such 3D cavities prove convincingly that the 
quality factor of low-contrast ring cavities can be raised 
enough to obtain devices with operation characteristics 
similar to those achieved in 2D designs. Theoretical and 
experimental work is ongoing on the optimization and 
demonstration of such an ultracompact integrated opti- 
cal circulator. 
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